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We analyze the diffractive (Pomeron) contribution to pp spin-fiip amplitude and discuss the pos- 
sible scenarios for energies available at the Relativistic Heavy-Ion Collider (RHIC). In particular, we 
show that RHIC data will be instrumental in assessing the real contribution of diffraction to spin 
amplitudes. 
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I. INTRODUCTION 



The study of diffraction on high-energy scattering using Regge formalism is a well-known subject Until the 
, coming in operation of Relativistic Heavy- Ion Collider (RHIC), the highest energy pp data available from accelerators 
were those of the Intersecting Storage Ring (ISR) at CERN. Those experiments, however, did not measure the 
polarization of the projectiles. Since the differential cross section [da/dt) at ISR energies is usually assumed to be 
dominated by the spin-non-flip amplitude, many of the models describing the diffraction in pp scattering do not take 
into account the contribution from spin-flip amplitudes |[ |^ . Others do (5[ 0, ||] but in this case they use data at 
^sj lower energies However, the study of diffraction in spin- flip amplitudes began many years before that when one of 
I us ]lO| noticed that polarization data at Tip scattering suggested a diffractive contribution in the spin-flip -np amplitude 
,. at high energy which becomes evident when the kinematical zero is removed. The spin-flip amplitude without the 
O ^' kinematical zero is named "reduced" and manifests itself as a peak in the forward direction which does not appear to 
' , vanish as the energy increases. In Ref. [ p^ it was explicitly noticed that, once the kinematical zero is removed, all 
partial waves act coherently in the small angle domain as it is typical of diffractive events . The following statement 
(-H ' was made later " the residual spin-flip amplitudes behave very much like spin-non-flip amplitudes at high energies 
• • I and exhibit a pronounced forward peak which is largely independent of the particular elastic reaction choseri^ . The 
same conclusion was obtained in the analysis of pp data some years later but in the case of pp scattering the 
situation is complicated by the existence of five independent helicity amplitudes, named 0j, with j = 1, . . . ,5. Of 
$H ' those, 01 and (j)^ are spin-non-flip amplitudes, <j)2 and (/14 are double spin-flip amplitudes, and ^5 is a single spin-flip 
amplitude. 

In this work we are interested in high energy and not too high t so that we can concentrate the analysis on the 
main aspects of the spin-flip amplitude for pp scattering in the diffractive region. In Ref. Q the magnitudes of 02 
and 04 with respect to the spin-non-flip amplitude were analyzed and reasonable arguments were given concerning 
the linear dependence on t for 02 (02 oc t when t — > 0) JlSf , while in the case of 04 the relation 04 oc i is consequence 
of angular momentum conservation. 

So, the general form of the polarization P for pp scattering 



[|0lP + |02P + |03P-H|04P-f4|05P] ^ ^ 



can be simplified assuming that 02 and 04 are small compared to the others amplitudes in the t region under interest. 
Using the definitions 

01 = — n — ' '?'5 = — 7, — ' (2) 
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where g{s,t) and h{s,t) will be considered effective spin-non-flip and spin-flip amplitudes, respectively, and assuming 
01 = (/"s the polarization can be rewritten as 

^ ^ Im[g{s,t)h*{s,t)] 
\g{s,t)\^ + 2\h{s,t)\^- 

Today, a rather considerable amount of data at higher energies has been gathered |^ in the relatively small angle 
domain and new perspectives are being opened by the coming in operation of the Relativistic Heavy Ion Collider 
(RHIC), the ideal machine to study polarization in high energy collision processes [p^ . 

In addition, our phenomenological information on the spin-non-flip amplitude is today much more complete and 
this can be used to reduce the uncertainties in the analysis. Using an explicit parametrization for the pp spin-non-flip 
amplitude fl^ whose parameters have been calculated against all high energy pp and pp data (except polarization) , 
we analyze the structure of the reduced spin-flip contribution, where the kinematical zero is removed with the factor 
\/^. The following conclusions are reached in this analysis. 

(a) The (reduced) spin-flip amplitude exhibits the typical peak in the forward direction which characterizes diffrac- 
tive amplitudes but the reduced spin-flip amplitude is less than 10% of the imaginary part of the spin-non-flip 
amplitude, (b) The data fitting with the same energy dependence in g{s, t) and /i(s, t) seems the best choice; a zero is 
automatically produced in the polarization; this moves with energy towards t = as a consequence of the analogous 
shift of the dip in da/dt induced by the zero of the spin- non-flip amplitude, (c) The magnitude of the polarization 
decreases as the energy increases but the extrapolation to 500 GeV predicts a non-negligible contribution if the same 
Pomeron trajectory for both spin-flip and spin-non-flip amplitudes is used. 



In the next section we present the expressions of the amplitudes and explain them. In Sec. [II we show the fit 
of the data and the modifications on the spin-flip amplitudes. We also consider a similar (nondiffractive) analysis of 
the data. Although the fit to the lower energy data is essentially the same, the two analysis predict a very different 
behavior at increasing energies. As we will stress, RHIC data should provide a reasonably clear cut answer to the 



question of whether or not diffraction contributes to the spin-flip amplitude. In Sec. IV we present the conclusions 



II. DEFINITION OF THE AMPLITUDES 

The effective pp spin-non-flip amplitude will be written as 

g{s,t)=a"f{s,t) = a+{s,t)-a^{s,t) (4) 

with 

a+{s,t) ^ ap{s,t) + af{s,t) and a^{s,t) — ao{s,t) + ai^{s,t), (5) 

where aip{s,t) and ao{s,t) are the Pomeron and Odderon amplitudes, respectively, and Of (s, t) [a^(s,t)] is the even 
(odd) secondary Reggeon These different amplitudes are taken directly from Ref. |l4j and their explicit forms 
are given in the Appendix together with the values of their parameters. 

To write the effective spin-flip amplitude h{s, t): (i) first, we neglect the spin contribution of secondary Reggeons to 
check whether diffraction could be the dominant phenomenon for the spin-flip amplitude at small |t| at the available 
energies; (ii) second, we notice that the combination of the exchange of two and three gluon ladders induce CP-even 
(and -odd) contributions to the spin-flip amplitude that affects Pomeron exchange. We have, therefore, to anticipate 
that the diffractive (Pomeron) contribution to the spin-ffip amplitude may have a slightly modifled structure from 
the usual one. Differently stated, to take into account three gluon ladders, we do not take a strict Regge pole 
parametrization for the Pomeron spin-flip term and we allow a real part contribution to be present. Thus, in the very 
small |t| domain, we write 



his,t) = a"-^(s,i) = (i7i+5i)^s"°'(*)e'3i*e(0.5-|f|) 

nip 

+ (z72 + <J2)^ S"°'(*^e''^*e(|i|-0.5), (6) 
nip 

where the mass of the proton rup is used to make the parameters dimensionless. In Eq. (H) 5 = (s/so)e^*'^/^, O is the 
step function and we assume sq = 1 GeV as in Ref. 0; the superscript sf (for "spin-flip") allows us to check if the 
IP trajectory can be different for spin-ffip and spin-non-ffip amplitudes. While the complex phase in Eq. (^) is used 
to take into account in a phenomenological way the mixing of CP even and odd contributions cited above, the step 
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function is used to account in the most economical way for the fact that the fitting procedure requires a change in the 
slope of h{s,t) (similar to what happens in da/dt) to have a good description of the polarization data. The precise 
point where the slope changes, however, is rather t insensitive and could be taken at almost any value between 0.2 
and 0.7 GeV^ (in We do not have an explanation for this effect nor for the insensitivity in the choice of the point 
where the slope changes but it proves that our parametrization is quite stable. To start with, we take the spin-flip 
Pomeron trajectory a^^ (t) to be exactly the same as derived for the spin-non-flip amplitude 

a'f (t) = ap (t) = ap (0) + a'pt, (7) 

where ap{Q) and a'p are given in the Appendix . 

III. FIT TO THE EXISTING pp POLARIZATION DATA 

The differential cross section da/dt for pp scattering is available at energies up to 63 GeV in the cm. system but 
polarization data are available only at lower energies. We utilized the data at = 13.8, 16.8, and 23.8 GeV (a 
total of 64 points) in the fit and we checked the quality of the result describing the polarization at 19.4 GeV. The 
parameters obtained in the fit are shown in Table | with the best P^r degree of freedom. In principle (and in 
practice) , the inclusion of the spin- flip amplitude in the game would require refitting all the parameters to reproduce 
the angular distributions as well. We will check, however, that the fit of the polarization data without performing a 
new fit to the angular distribution does not modify the quality of the fit to the latter while making much more direct 
the analysis of the polarization data 

Figure ^ presents the set of polarization data used in the fit together with our reconstruction. Figure |^ shows the 
polarization at ^/s = 19.4 GeV with the prediction of our model since this set was not used in the fit. In Fig. ^ 
we show da/dt at various energies as described by adding the (squared) spin- flip amplitude [see Eq. ( |A.9D in the 
Appendix ] to the (squared) spin non flip term. As anticipated, the quality of the fit has not been altered. 

Several comments are in order. 

(a) The IP contribution to the spin-flip amplitude h{s, t) is considerably smaller than to the spin-non-flip term g(s, t) 
(about 5%, roughly, at ^/s = 20 GeV and changing little until 500 GeV). This value is compatible with the analysis 
performed in Q] where use was made of the relative amplitude — to05/(\/— i lm(j)^). In our case, Im rs = —0.054 
at ^/s = 500 GeV. We should mention that we find natural that the spin-flip part of the amplitude should be a small 
fraction of the spin-non-flip but, again, we have no real reason for this, (b) Contrary to the discussion made in Ref. 
| p2t , the (small \t\) slope of the spin- flip amplitude /3i = 4.74 GeV~^ appears to be not very much different from the 
effective slope of the Pomeron in spin- non- flip amplitude (see the Appendix ). The present parametrization, however, 
is considerably more elaborate and the set of data corresponds to higher energies so the results of these papers may 
not be directly comparable, (c) Our result [i.e., h{s,t)] cannot be extended to \t\ values much higher than few GeV^ 
because the spin-non-flip amplitude utilized is valid at the Born level [ p^ and its description in the region after the 
dip {\t\ > 1.5 GeV^) is not very good. To extend our considerations to higher \t\, it would be necessary to adopt the 
more sophisticated eikonalized version. Anyway, the t region of interest for RHIC is up to 1.5 GeV^ H so we can 
confine our analysis to the not too high \t\ region. This point may have to be reconsidered in the future along with 
more detailed analysis, (d) The spin- non- flip amplitude in Ref. [Q fits da/dt without the spin contribution and that 
description is not spoiled by the presence of h{s, t) on this work since it is very smaller than g(s, t). 

A. The kinematical zero of the spin-flip amplitude 

The original work about diffraction and polarization data at -np scattering calculated the reduced spin-flip amplitude 
removing the kinematical zero by means of a sin 6 instead of [|l^ . Since the energies available at that time were 
much lower, it would be difficult to see differences between the two approaches. But the data presently available and 
the coming in operation of RHIC at even higher energies raises the question of what would happen if the factor sin 9 
was used. In principle, the use of relativistically invariant variables used earlier seems more appropriate but we feel 
that the answer can only come from experiments. 

To answer this question we adopt in this section the spin-flip amplitude 



h{s,t) = a^^(s,t) = (i7i +(5i)sin6l s" (*)e''i*e(0.5- |t|) 
+ (i72 + (52)sin6'I"''^(*)e''=*e(|i| -0.5). 



(8) 



4 



We use the same procedure of the fit with Eq. that is, we fit the data at ^/s = 13.8, 16.8 and 23.8 GeV assuming 
a^^ {t) = a]p{t). Then we check the quaUty of the fit with the description of the polarization data at 19.4 GeV. The 
values of the parameters obtained with Eq. (||) are shown in Table ||. 

Following the same procedure of the previous section, we present in Fig. ^ the result of the fitting while Fig. ^ 
shows the polarization at ^/s = 19.4 GeV and Fig. || presents da/dt. Two differences in the present results and 
in those of the previous analysis are worth being emphasized. First, the ratio of the reduced spin-flip [h(s , t) / sin 6] 
to the spin-non-fiip amplitude is much bigger (above 90%) and this big ratio was noticed at lower energies too [ p^ . 
Second, the slope f3i — 6.25 GeV^ is higher than the previous result (Table ||) so the reduced spin-flip amplitude is 
steeper on this case. 

Although Eqs. (^ and (||) present very different values for the parameters Si, and 71, Figs. |l|-|| show that 
the descriptions of the data are similar and the fits have a comparable x^- To understand the consequences of using 
one parametrization or the other it is necessary to remember that siu.0 oc \/ —tj s. So, to all effects, Eq. (^) does 
not have the s dependence of the Pomeron. Differently stated, the parametrization (P), ultimately, is not diffractive. 
The comparison between the two parametrizations used for h{s,t), therefore, Eqs. and (^ goes to the root of 
the problem of whether or not diffraction contributed to spin-flip. This difference will have relevant consequences at 
higher energies and the pp2pp experiment at RHIC JTsI will be able to discriminate between these two scenarios. 

B. Predictions to higher energies 

The pp2pp experiment at RHIC will provide information about polarization in pp scattering at energies between 50 
and 500 GeV 13|. That means a big increase in the amount of information about the spin content of the proton so 
it is important to see what we can say about the polarization in that energy range based on the information we have 
derived from the data available at lower energies. To this aim, we calculate the polarization at ^/s = 50 and 500 GeV 
with Eq. (^ plotting the result in Fig. ^. We can see that the polarization decreases in magnitude with increasing 
energy but it is still sizeable at 500 GeV. Equation (^) produces a peak of about 10% of positive polarization around 
—t — 1.25 GeV^ and ^/s = 500 GeV as compared with a (positive) polarization of over 15% at y/s = 50 GeV. 

Now we compare this result with the prediction using Eq. (||) (see Fig. H). The polarization calculated with the 
factor sin 9 is much smaller (about 5% at y/s = 50 GcV) but becomes essentially zero at y/s = 500 GeV. As already 
mentioned the factor l/y/s hidden inside the sine function is important to separate the predictions at high energies. 

We have also checked the possibility of a different intercept for the trajectory a^^ (t). In this case Eq. (0) is modifled 
and a^-^ {0) becomes a new parameter to be fltted together with the other six parameters. The values obtained (with 
a similar /^df) for them were quite absurd, including a negative intercept a''''^(0) when the factor y/^/nip is used, 
and the polarization predictions at RHIC energies are so small that it would be essentially impossible to detect those 
values experimentally. Since there is no improvement of the flt with the increased number of parameters we discard 
this solution. 



IV. CONCLUSIONS 

We analyzed the spin-flip amplitude and removed its kinematical zero to study the Pomeron contribution to the 
spin as suggested long time ago |T^. We made use of two hypothesis to remove the zero, ^/^/rnp and sinfl, where 
the first is more convenient at high energies since it is relativistic invariant while the second was used in the original 
work about diffraction in pp spin- flip amplitudes [ p^ . The differences resulting from the application of Eqs. (^ or 
(^ are quite large although the descriptions of the data between ^/s — 13 and 24 GeV are very similar. The fraction 
of h{s, t) to g(s, t) is small if \/^/mp is used (around 5%) but it becomes big when sinO is utilized (more than 90%). 
Although, as already mentioned, we do not have a strict reason to prefer one over the other, the flrst solution is much 
more in line with the traditional analyses. Also, the slope at small \t\ changes when ^/^/rrip is substituted by sin 6*, 
showing that h{s,t) becomes steeper with the second option. At the same time, the slope f32 is practically the same 
in both cases showing that there is a modification on h{s,t) around —t = 0.5 GeV^, which conflrms the necessity of 
two slopes in the spin-flip amplitude as it was assumed. The extrapolations shown in Figs. |^ and ^ suggest that the 
data from RHIC will be crucial to understand and choose the best form to describe the spin-flip amplitude. 

In summary, the prediction is very straightforward: a positive polarization of the order of 10% is predicted at RHIC 
energies of ^/s — 500 GeV if the Pomeron (diffraction) contributes to the spin-flip amplitude as lower energy data 
tend to suggest. A vanishing contribution is otherwise expected. 

It would be possible to improve the description of the data with more elaborated forms for both spin-non-flip and 
spin-flip amplitudes (by eikonalizing them, for example, Ref. [ p4| , by introducing subasymptotic effects in the spin-flip 
part by secondary Reggeons and so on) but we expect that the main conclusions of our analysis would remain the same 
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in the region of small t (where the Born amplitudes work well) and high energies (where the Pomeron dominates). A 
more detailed analysis to separate the contributions of the exchanges of two and three gluon ladders would be more 
gratifying from the theoretical point of view as well as a better study of the break in the diffractive slopes. All these 
points will, hopefully, be reconsidered in the future. 
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APPENDIX: THE SPIN-NON-FLIP AMPLITUDE 

The spin-non-fiip amplitude utilized in this work is 

a'^^{s,t) = a+{s,t) - a_{s,t), (A.l) 

where 

a+{s,t) — a]p{s,t) + af{s,t) (A. 2) 

and 

a-{s,t) — ao{s,t) + au){s,t). (A. 3) 

The expressions for the two Reggeons used in Rcf. are 

aii(s,i) =aflS"«^*^e^«* (A.4) 

and 

anit) = aniO) + a'j^t {R ^ f and w) (A.5) 

with af{auj) real (imaginary). 

For the Pomeron, the spin-non-flip amplitude is 

(s, t) = aj.r^(*) [g6^[a^(t)-i] + In 5) + dp In s] (A.6) 

while for the Odderon we use 

ao(s,t) - [l-exp(7t)]ao5"°(*)[e''°("°(*)-i)(&o+ln.5) + dolnS], (A.7) 

where again a]p{ao) real (imaginary) and we have utilized ai{t) = ai(0) -I- a^<, where i = IP^O. 
Our deflnition for the amplitude follows H| so that 



at = — Im{a"^(s,i = 0)}, (A., 
s 



da 



=^[\a-f{s,t)\^ + 2Wf{s,t)\\ (A.9) 

In this work we retain the same parameters for the spin-non-flip amplitude as in Ref. Jl^ and we keep them fixed 
while fitting the parameters of the spin-flip amplitude. We utilize the dipole model at the Born level since most of 
the polarization data is contained in the t domain corresponding to the region before the dip in da/dt (well described 



without eikonalization). The values of the parameters of the spin-non-fiip amplitude |14| are shown in Table III. 
To calculate the polarization we utilized the form 

lm{a-fis,t)[a^f{s,t)r} 

|a"/(s,t)|2-f 2|a^/(s,i)|2' ^ 
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where the star on the numerator means complex conjugate. 
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TABLE I: Values of the parameters obtained from fitting polarization data at = 13.8, 16.8, and 23.8 GeV with Eqs. (|^) 
and (0). The xV^df is 1.1. 



i 


7i Si 


A (GeV-") 


1 


1.35 X 10"^ 2.64 X 10"^ 


4.74 


2 


2.55 X 10"^ 5.38 x 10"^ 


2.29 


TABLE II: Resuhs from fitting polarization data at = 13.8, 16.8, and 23.8 GeV with Eqs. (|) and (||). 


The xV^DF is 1.1. 


i 


7i St 


Pi (GeV-") 


1 


2.55 4.80 


6.25 


2 


0.18 0.45 


2.30 


TABLE III: Parameters of the dipole model at the Born level [Q with i = ]P,0, f, uj. 




Pomeron Odderon / Reggeon 


u) Reggeon 


«40) 


1.071 1.0 0.72 


0.46 


a'. (GeV- 2) 


0.28 0.12 0.50 


0.50 


Oi 


-0.066 0.100 -14.0 


9.0 


h 


14.56 28.10 1.64 GeV"^ 


0.38 GeV"" 


di 


0.07 -0.06 




7 (GeV"^) 


1.56 





0.10 




-0.30 - 





FIG. 1: Results from fitting polarization data at (a) 13.8 GeV, (b) 16.8 GeV, and (c) 23.8 GeV (see Table 
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FIG. 2: The prediction for the polarization at 19.4 GeV (not used in the fit) compared with the experimental data at that 
energy (parameters from Table |). 
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FIG. 3: The differential cross section obtained in this work taking into account the spin-flip amplitude, Eq. (fj). The highest 
set of data corresponds to 23.5 and 27.4 GeV grouped together. The other sets (multiplied by powers of 10~^ j are 30.5, 44.6, 
52.8, and 62 GeV. 




FIG. 4: Results from fitting polarization data at (a) 13.8 GeV, (b) 16.8 GeV, and (c) 23.8 GeV (see Table 
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FIG. 5: The prediction for the polarization at 19.4 GeV (not used in the fit) compared with the experimental data at that 
energy (parameters from Table nWj. 
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FIG. 6: The differential cross section obtained in this work taking into account the spin- flip amplitude through Eq. (g). About 
the data see Fig. H. 
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FIG. 7: The polarization predictions for — 50 and 500 GeV with the parameters of Table |. 
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olarization predictions for 50 and 500 GeV with a detailed view of the data at 500 GeV in the inset (parameters 



